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ABSTRACT 

The  studies  covered  in  this  report  are  the  design 
of  optical  frequency  shifters  for  AC  interferometry 
(Section  II)  and  the  production  and  assessment  of  corrector 
plates  used  to  optically  compensate  for  distortions  in 
laser  beams  (Section  III) . 

Section  II  describes  efforts  to  work  out  an  electro- 
optical  system  for  producing  Doppler  frequency  shifts  in 
laser  beams.  Earlier  efforts  to  build  practical  AC  inter¬ 
ferometers  ran  into  difficulties  with  the  frequency 
shifters  required:  vibrating  mirrors  and  other  mechanically 
moving  frequency  shifters  had  to  be  discarded.  Early  ex¬ 
periments  with  electro-optical  frequency  shifters  had 
practical  difficulties  in  obtaining  high  enough  spectral 
purity  in  the  output  beam.  These  difficulties  have  now 
been  overcome. 

Section  III  describes  the  design,  fabrication,  and 
testing  of  a  correction  plate  for  removing  distortions 
in  a  laser  beam.  The  work  consisted  of  1)  measuring  the 
deformation  to  be  corrected,  2)  digitizing  the  deformation 
information,  3)  preparing  a  paper  tape  to  control  an 
optical  milling  machine,  4)  machining  a  scaled-up  correction 
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plate,  and  5)  using  the  correction  plate  with  the 
appropriate  index-matching  fluid  to  attempt  to  correct 
the  original  deformation.  The  pilot  experiment  clearly 
demonstrated  the  feasibility  of  this  approach  to  correct¬ 
ing  arbitrary  wavefront  deformations.  With  our  current 
capability  we  could  easily  reduce  wavefront  distortions 
of  as  much  as  50  wavelengths  to  less  than  one  wave¬ 
length.  Our  approach  would  be  particularly  well-suited 
to  the  correction  of  large  numbers  of  optical  components, 
since  all  required  operations  can  be  done  relatively  in¬ 
expensively  and  routinely  after  the  initial  set-up. 
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I .  REPORT  SUMMARY 


The  work  described  in  this  report  is  a  continuation 
of  earlier  studies  under  Contract  F33615-70-C-1451  and 
covers  two  principal  areas  of  interest;  design  of  an 
electro-optical  frequency  shifter  for  use  in  AC  interfero¬ 
metry,  and  the  production  and  testing  of  a  scaled  correction 
plate  for  optical  correction  of  beam  divergence  in  crystal 
and  glass  laser  rods. 

Design  of  an  Electro-Optical  Frequency  Shifter 

We  found  during  the  course  of  our  experimental  in¬ 
vestigation  that  a  primary  problem  in  building  practical 
AC  interferometers  was  obtaining  frequency  shifters  that 
could  be  used  to  produce  tbs  local  oscillator  beams  needed 
for  them.  Most  of  our  early  experiments  utilized  vibrating 
mirrors  to  produce  frequency  shifts  in  laser  beams  by  Doppler 
methods  (a  technique  called  serrodyning) .  We  found,  however, 
that  it  was  difficult  to  obtain  reliable  operation  from  such 
devices,  and  that  more  sophisticated  frequency  shifters  were 
needed.  At  the  same  time,  the  bandwidth  requirements  of  our 
system  prohibited  the  use  of  mechanically  moving  frequency 
shifters,  so  we  resorted  to  an  electro-optical  system. 

In  our  early  experiments  with  an  electro-optical  frequency 
shifter,  we  found  that  the  spectral  purity  of  the  output  beam 
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from  the  frequency  shifter  was  not  sufficiently  high  to 
make  it  useful  for  an  AC  interferometer.  We  therefore  set 
about  a  series  of  experiments  to  determine  the  effect  of 
various  modulator  characteristics  on  the  spectral  purity 
of  a  frequency  shifted  output  beam.  Both  single-pass  and 
double-pass  configurations  were  tried.  The  output  beam 
had  greater  spectral  purity  in  the  double-pass  configuration. 

AJlethod  of  Optically  Correcting  Crystal  and  Glass  Laser  Rods 

The  beam  distortion  to  be  corrected  is  determined  by 
obtaining  an  interferogram  of  a  wavefront  which  has  been 
distorted  by  passage  through  the  laser  rod  responsible  for 
the  distortion.  This  interferogram  is  essentially  a  contour 
map  of  the  distorted  wavefront,  and  also  a  contour  map  of 
the  desired  optical  thickness  of  the  finished  correction 
plate.  The  plate  is  made  to  conform  to  the  interferogram 
contours  except  that  high  spots  on  the  interferogram  become 
low  spots  on  the  correction  plate,  and  vice  versa.  If  the 
original  beam  is  passed  through  this  correction  plate,  a 
distortion-free  beam  should  result. 

For  accuracy  in  machine  contouring  of  the  plate  surface, 
we  enlarge  the  distortion  by  a  scale  factor  M  .  Contour¬ 
ing  is  done  semi -automatically  by  an  optical  milling  machine 
controlled  in  depth-of-cut  by  a  paper  tape  generated  using 
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data  taken  from  a  photograph  of  the  interf erogram .  The 
contoured  ulate  must  then  be  scaled  back  down  to  the  proper 
magnitude  tor  use.  This  is  done  by  immersing  it  in  a  fluid 
of  refractive  index  nearly  matching  the  index  of  the  plate. 

An  index  mismatch  of  An  produces  a  scale  factor  M  =  1/  r. 
The  near -matching  technique  also  does  away  with  the  necessity 
for  optically  polishing  the  correction  plate,  since  the 
surface  roughness  is  scaled  down  by  the  same  factor  M 

The  final  steps  in  the  study  involve  immersion  of  the 
completed  correction  plate  in  the  near-matching  index  fluid 
and  testing  of  the  original  beam  for  absence  of  distortion. 

The  correction  plate  proved  to  be  successful.  When  the  plate 
was  inserted  into  the  beam,  a  marked  improvement  in  the 
quality  of  the  transmitted  wavefront  was  obtained:  20  fringes 
of  deformation  in  the  uncorrected  beam  were  reduced  to  only 
a  few  fringes.  Even  better  results  can  be  achieved  with  a 
more  refined  version  of  the  same  basic  technique. 


II.  DESIGN  OF  OPTICAL  FREQUENCY  SHIFTERS 
FOR  AC  INTERFEROMETRY 

1 .  Single  Sideband  Suppressed  Carrier  Modulation  of  Light 
by  Electro-Optic  Techniques 

The  method  of  detecting  a  signal  by  mixing  it  with 
a  reference  signal  of  different  frequency  (heterodyning) 
has  been  widely  used  in  the  domain  of  radio  communications. 
Since  the  advent  of  lasers,  the  properties  of  heterodyne 
detection  have  been  applied  to  optical  communications  as 
well,  and  seem  to  have  promising  applications  in  such 
areas  as  Laser  Doppler  Velocimetry  and  fringe-counting 
interferometry,  where  the  heterodyne  signal  is  easily 
related  to  the  speed  of  the  moving  particles  or  the  dis¬ 
placement  of  an  object  under  test. 

The  Doppler  shifting  of  optical  frequencies  can  be 
produced  by  several  means.  Most  of  those  now  in  use  are 
mechanical;  for  example,  the  linear  translation  of  a 
mirror,  linear  rotation  or  translation  of  a  grating,  ro¬ 
tation  of  a  half  wave  plate,  etc.  Other  methods  are 
electromagnetic  in  nature,  using  the  change  in  refractive 
index  of  the  medium  due  to  magnetic,  acoustic  or  electric 
fields.  Although  the  mechanical  means  have  the  advantage 
of  ease  in  alignment  and  simplicity,  their  use  is  limited 
by  their  very  low  bandwidth,  especially  in  interfe  ometric 
applications  where  they  can  upset  the  experiment  by  the 
vibrations  they  generate.  Magneto  optic  methods  (Zeeman 
effect)  require  rather  high  magnetic  fields,  because  of 

Preceding  page  blank 
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the  smallness  of  the  Verdet  constant  for  most  media. 
Acousto-optic  frequency  shifting  (Bragg  diffraction) 
usually  alters  the  direction  of  the  beam,  which  makes 
it  very  useful  for  constructing  beam  deflectors  or  holo¬ 
graphic  memory  readers,  but  unfit  for  the  interferometric 
applications  we  are  mainly  interested  in. 

Hence,  we  will  limit  our  attention  to  electro¬ 
optic  frequency  shifters,  comparing  them  in  efficiency, 
power  requirements,  ease  of  adjustment,  etc. 


' 

2.  Theory  of  Electro-Optic  Frequency  Shifting 

The  electro-oot ic  modulation  of  light  makes  use  of 
the  fact  that  two  different  polarizations  of  light  exper¬ 
ience  two  slightly  different  indices  of  refraction  when 
going  through  an  electro-optic  crystal  upon  which  some 
voltage  is  applied.  By  modulating  this  voltage,  one 
changes  the  optical  path  through  the  crystal  at  will,  and, 
hence,  one  modifies  the  instantaneous  frequency  of  the 
output  light. 


Vo//*aje  V*V(0 


The  frequency  of  the  output  light  depends  on  the 
modulating  function  V(t),  and  in  the  following  pages  we 
will  investigate  several  of  the  voltage  excitations  that 
must  be  used  in  order  to  obtain  a  constant  shift  in  the 
frequency  of  the  input  light. 

Without  going  too  far  back  into  the  general  prop¬ 
erties  of  electro-optic  crystals,  we  recall  that  the 
dielectric  constant  of  a  homogeneous,  nonconducting  and 
anisotropic  medium  is  a  symmetric  tensor,  so  that  it  is 
made  up  of  6  independent  components.  The  structure  of 


an  anisotropic  medium  allows  two  monochromatic  plane  waves 
with  two  different  polarizations  and  two  different  veloc¬ 
ities  to  propagate  in  any  given  direction,  and  both  the 
velocities  and  the  directions  of  polarization  can  be  easily 
found  from  the  index  ellipsoid  of  the  medium  built  from 
the  dielectric  tensor: 
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where  the  standard  notation  is  used: 
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The  two  planes  of  polarization  which  are  allowed 
in  the  crystal  to  an  incoming  beam  are  obtained  by  consi¬ 
dering  the  intersection  of  the  index  ellipsoid  with  the 
plane  perpendicular  to  the  direction  of  propagation  of 
the  beam:  if  it  is  an  ellipse,  the  allowed  polarizations 
are  parallel  to  the  main  axes  of  the  ellipse  and  the 
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corresponding  indices  of  refraction  are  half  the  length 
of  these  axes.  If  the  intersection  is  a  circle,  then  the 
crystal  is  isotropic  for  that  particular  direction  of 
propagation.  When  the  crystal  is  subjected  to  an  electric 
field,  the  die1ectric  constant,  and  hence,  the  index  ellip¬ 


soid  undergo  a  change  such  that: 


''  'wM 
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I  /w» 
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where  the  Ej’s  are  the  three  components  of  the  field  on 
the  x,  y,  z  axes. 

Usually,  for  most  modulation  purposes,  the  crystal 
is  cut  in  such  a  way  that  when  no  voltage  is  applied,  the 
crystal  is  isotropic  for  the  incoming  light.  The  setting 
of  an  electric  field  makes  the  crystal  anisotrcDic,  the 
directions  of  the  fast  and  slow  axes  depending  on  the  di¬ 
rection  of  the  electric  field  and  the  symmetry  of  the 
crystal.  The  retardation  introduced  between  the  two 
polarizations  allowed  to  propagate  through  the  crystal 
is  simply: 


r  (VO 


where  d  is  the  thickness  of  the  crystal.  The  crystal  acts 
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then  as  an  r-wave  retarder,  where  the  retardation  r  is 
proportional  to  the  applied  field  (linear  Pockels  effect) . 


A.  Case  of  One  Crystal 
1.  Transverse  field. 

The  effect  of  a  rotating  r-wave  plate  upon  a  circu¬ 
larly  polarized  beam  of  light  is  described  by  the  Jones 
Calculus.  The  matrix  corresponding  to  an  r-wave  plate, 
whose  axer;  are  at  angle  a  with  the  x-y  axes  is  the  follow¬ 
ing  : 


T 

r-  wave. 


cos2* 
sm  <x 


COS«<Sltto< 
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infi\  . 

Jsinx  cos* 
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The  action  of  such  a  plate  on  a  circularly  polarized 


beam  of  ligh-.  results  in  the  following  output  beam: 


Hence,  if  we  set  the  r-wave  plate  in  a  uniform  rota¬ 
tion,  the  output  light  is  made  up  of  two  components:  one, 
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with  amplitude  cosTrr  is  of  the  same  frequency  and  state 
of  polarization  as  the  input  light,  and  the  other  one, 
with  amplitude  sinfrr,  is  shifted  in  frequency  by  twice 
the  frequency  of  rotation  of  the  plate,  and  has  polari¬ 
zation  orthogonal  to  that  of  the  input  light.  If  the 
plate  in  use  is  a  half-wave  plate,  a  left  circular  beam 
of  light  is  totally  transformed  into  a  right  circularly 
polarized  beam,  up  or  down  shifted  in  frequency  (depen¬ 
ding  upon  whether  the  rotation  of  the  plate  is  in  the 
same  or  opposite  direction  as  the  direction  of  polari¬ 
zation  of  the  incoming  „ ight) . 

This  method  of  a  rotating  half-wave  plate  has  been 
used  for  phase  measurements  in  interferometry1-  and  micros¬ 
copy2.  However,  its  range  is  limited  to  frequencies  up 
to  several  hundred  Hertz,  the  mechanical  rotation  may 
introduce  vibrations  or  turbulence  in  the  interferometric 
set-up,  and  care  must  be  taker  \_o  compensate  for  a 
possible  uneven  thickness  of  the  plate,  which  would  gen¬ 
erate  a  supplementary  phase  modulation  of  the  beam. 

It  has  been  shown  theoretically  and  experimentally^ »  4 , 5 
that  the  same  effect  can  be  obtained  without  any  mechan¬ 
ical  rotation,  by  applying  a  rotating  transverse  electric 
field  to  crystals  of  cubic  symmetry,  and  of  point  group 
43m,  ^hen  light  is  being  propagated  along  the  threefold 
axis.  When  no  voltage  is  applied,  the  crystal  is  isotro¬ 
pic,  and  when  a  transverse  field  is  on,  the  amount  of 
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birefringence  induced  in  the  crystal  depends  onl>  upon 
the  magnitude  of  the  field,  but  not  upon  its  direction, 
although  the  direction  of  the  axes  of  birefringence 
are  related  to  the  direction  of  the  field,  as  it  can  be 

shown  in  the  next  lines,  taken  from  Kamikow  and  Turner's 
article"* : 

From  the  symmetry  of  a  cubic,  43m  crystal,  it  follows 
that  the  electro-optic  tensor  of  such  a  crystal  (for  exam¬ 
ple,  zinc  sulfide,  or  lithium  niobate)  is  such  that: 


When  the  transverse  electric  field  E  (E  ,E  )  is 

-  x  y 

applied  perpendicularly  to  the  direction  of  propagation 
of  light,  the  equation  of  the  index  ellipsoid  in  the 
z  =  0  plane  becomes: 


-4- 


this  ellipse  becomes; 
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\.here  E  is  the  amplitude  of  the  electric  field: 

Ex  =  E  cos  ^ 

F  =  E  si n  0 

7 

The  angle  between  the  principal  axes  of  the  ellipse 
and  the  x-y  axes  is: 


Hence,  as  the  transverse  electric  field  is  set  into 
rotation,  the  axes  of  induced  birefringence  rotate  twice 
as  slowly,  and  the  shift  in  the  frequency  of  the  optical 
beam  is  equal  to  the  frequency  of  rotation  of  the  electric 
field . 

Experimentally,  the  transverse  rotating  field  is 
obtained  by  feeding  two  pairs  of  orthogonal  plane  elec¬ 
trodes  respectively  with  sine  and  cosine  wave  voltages. 

The  method  has  been  tried  in  the  audio  range  with  a  zinc 
sulfide  crystal,  unfortunately  of  very  poor  optical  qual¬ 
ity^,  and  with  better  results,  at  110  MHz  with  li_hium 
niobate^ . 
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2.  Transverse  or  Longitudinal  Excitation. 

The  effect  equivalent  to  the  linear  translation  of 
a  mirror  to  shift  the  frequency  of  the  incoming  beam  may 
be  obtained  by  applying  a  linear  voltage  to  electro-optic 
crystals.  This  effect  (called  serrodyning)  has  been  used 
to  shift  the  frequency  of  gaussian  pulses^  by  synchro¬ 
nizing  laser  pulses  with  the  zero  crossings  of  an  applied 
sinusoidal  field.  It  can  also  be  shown^  that,  wr.en  this 
method  is  used  with  a  quasimonochromatic  input  light,  the 
total  energy  of  the  output  beam  is  found  at  the  shitted 
frequency  if  the  saw  tooth  modulation  is  used  and  the 
maximum  voltage  of  the  ramp  is  the  half-wave  voltage  of 
the  crystal. 

B.  Case  of  Two  Crystals 

Instead  of  using  the  transverse  electro-optic 
effect  or  linear  excitation,  as  in  serrodyning,  we  in¬ 
vestigate  the  possibilities  of  a  longitudinal  sine-wave 
excitation,  which  allows  the  use  of  good  optical  quality 
crystals  of  the  KDP  type.  These  crystals  are  quite  strain 
free,  and  do  not  become  opaque  after  being  exposed  to  U.V. 
radiation  as  LiNbC>3  does.  Because  they  absorb  very  little 
light,  the  problem  of  heating  is  avoided.  Since  the  ex¬ 
citation  is  along  the  direction  of  propagation  of  light, 
the  principal  axes  of  induced  birefringence  are  no  longer 
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determined  by  the  direction  of  the  electric  field,  but 
by  the  symmetry  of  the  crystal  itself.  Since  no  physical 
rotation  of  the  crystal  is  desirable,  one  has  to  use  more 
than  one  crystal  along  the  path  of  light  in  oroer  to 
obtain  the  frequency  shift  of  the  incident  light. 

Consider  rhe  case  of  two  Z-cut  crystals  in  the  longi¬ 
tudinal  mode,  so  that  their  optic  axis  is  parallel  to  the 
direction  of  propagation  of  the  incident  beam  and  they 
are  isotropic  when  no  voltage  is  applied.  Their  princi¬ 
pal  axes  are  at  an  angle  a  with  each  other.  The  voltages 
are  sine  wave  functions,  and  we  suppose  that  the  incoming 
light  is  elliptically  polarized: 


V,  a  Ve  cos,  t 

X=  V0  cos(uJ+y) 


And  = 


When  a  voltage  is  applied  longitudinally  on  a 
crystal,  the  index  ellipsoid  is  changed  from: 


Z ,  2. 

X  +y  +  2 


n2- 


to : 


n 


r-  +  -nr+27«E*+-fr  =  1 

n,  n. 
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The  changes  in  index  of  refraction  along  the  induced  fast 
and  slow  axes  are  then: 


O 


and  the  retardation  r  between  the  two  components  of  an 
optical  field  along  th^ se  axes  is: 


Since 


$  =  lizr  = 

An  =  i*  A 


u 


the  retardation  is  then: 

— -  nlr„  E.Js2JLn\V 


(,b  2 


6b 


Knowing  the  characteristics  of  the  input  light,  we 
deduce  easily  the  components  of  the  output  light  by  use 
of  the  transfer  matrices  of  the  two  crystals: 
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where  6t,  and  are  the  maximum  retardations  induced 
respectively  in  the  first  and  second  crystals,  and  are 
directly  proportional  to  the  maximum  voltages  applied 
to  the  crystals,  as  seen  from  the  relation  between  6 
and  V,  above. 

From  the  general  expression  giving  the  components 
of  the  output  light,  as  stated  above,  we  determine  which 
values  of  the  parameters  (i.e.  voltage  phase  shift,  angle 
between  the  principal  axes  of  the  two  crystals,  voltage 
amplitude)  maximize  the  output  energy  in  the  desired 
sideband.  We  consider  that  the  incident  light  is  circu¬ 
larly  polarized,  and  that  the  voltages  applied  to  both 


crystals  have  the  same  amplitude.  Then,  the  output  can 
be  written  air  ^  r 


E  =  £ 


xren  a:  r 

t  +  CO  S  -iA.Jcosa,f--cos(oJ^)/ 


COS*  -  L  £ 


sS  tn<x 


COS  tom 


far  ^f)]  '  -  i  ^os  />  cos[(jn  hf)J 


Sir}  c*  -  L  £ 


COS* 
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The  spectral  content  of  the  outp;  t  depends  only 
upon  6 , <J>  and  a,  and  the  following  graphs  describe  the 
behavior  of  the  first  harmonics  in  the  output  for.  in¬ 
creasing  values  of  the  induced  retardation  6,  with  respect 
to  the  two  other  parameters.  An  analysis  of  Figs.  1  and  2 
shows  that  complete  elimination  of  the  lower  sideband  is 
obtained  when  the  electrical  voltages  applied  to  the 
crystals  are  in  quadrature  (4>  =  and  the  output  of  the 
desired  sideband  is  maximized  when  the  principal  axes 
of  the  crystals  are  at  45°  with  respect  to  each  other. 

These  conclusions  are  true whatever  the  maximum  induced 
retardation  is. 

On*  may  note  that  these  results  are  reversedr  that 
is,  the  upper  sideband  is  totally  eliminated  and  the  out¬ 
put  of  the  lower  sideband  is  maximized  if  the  two  crystals 
are  rotated  by  instead  of  -J-  ,  or  if  the  phase  difference 

between  the  two  voltages  has  its  sign  changed.  The  lower 
sideband  can  also  be  maximized  if  right  circularly  pol¬ 
arized  light  (  instead  of  left  circularly  polarized)  is 
incident  upon  the  original  set-up.  Hence,  one  may  elim¬ 
inate  at  will  one  of  the  two  sidebands  and  maximize  the 
other  one  by  one  of  the  following  changes: 

-  right  circular  polarization  - ►  left  circular 
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tilde  3**  harmonica 


Spectral  content  of  the  output  versus  the  phase  shift  $  between  the  voltages  applied 


Let  the  incident  light  be  circularly  polarized,  the 


voltages  applied  to  the  crystals  be  of  equal  amplitude 
and  in  quadrature,  and  the  crystals  oriented  in  such  a 
way  that  their  birefringent  axes  are  at  45  degrees  with 
respect  to  each  other.  Then  the  Jones  vector  correspon¬ 
ding  to  the  output  light  is: 


Using  the  following  relations: 

izco%6  y  k  /  \  ik8 
£.  =  Zj  i 

Tk(-z)  =  f)k  Jt(z) 


one  finds  that  the  components  of  the  output  light  are: 


The  output  beam  can  then  be  decomposed  into  an 
infinite  series  of  right  and  left  circularly  polarized 
harmonics,  the  even  ones  having  the  same  state  of  polari¬ 
zation  as  the  incident  beam,  the  odd  ones  being  orthogonal 
to  it,  and  with  intensities  given  by: 


By  letting  the  beam  go  through  a  circular  polarizer, 
one  eliminates  all  the  even  harmonics,  and,  since  on  electro¬ 
optic  crystal  is  always  operated  with  voltages  no  higher 
than  its  half-wave  voltage,  the  first  harmonic  is  always 
much  larger  than  the  other  odd  harmonics  exiting  from  the 
circular  polarizer,  and  to  a  first  approximation,  the  inci¬ 
dent  energy  of  the  beam  at  frequency  too  is  transmitted 
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through  the  system  after  being  shifted  to  the  frequency 


Experimentally,  we  have  set  up  this  type  of  electro¬ 
optic  frequency  shifter  by  using  two  potassium  dideuteriu 
phosphate  cells  (KD*P) ,  whose  half-wave  voltage  is  about 
1800V,  because  each  is  made  of  two  crystals  optically  in 
series  and  electrically  in  parallel.  In  order  to  compare 
the  output  of  the  system  to  the  input  light,  we  insert  it 
in  a  Mach  Zehnder  interferometer: 


The  linearly  polarized  light  from  the  HeNe  laser  is 
made  circularly  polarized  by  means  of  the  first  quarter 
wave  plate  (QWP1).  it  then  goes  through  the  two  electro¬ 
optic  cells,  Cl  and  C2,  and  the  circular  polarizer,  made 
up  of  a  quarter  wave  plate  followed  by  a  linear  polarizer 
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This  final  assembly  brings  complete  extinction  when  no 
voltage  is  fed  to  the  cells.  This  voltage  is  generated 
by  meand  of  an  oscillator,  a  phase  shifter  (PS)  with  two 
outputs  whose  relative  phases  can  be  adjusted,  a  20  W 
stereo  amplifier  and  two  step-up  transform, ers  (Tl  and  T2)  . 

In  order  to  evaluate  the  performance  of  the  actual 
set-up,  we  derive  the  output  spectrum  of  the  frequency 
shifter  when  some  of  the  components  are  misaligned  or 
present  some  defects,  or  when  the  input  light  is  not 
circularly  polarized. 

For  example,  we  show  the  computations  involved  in 
two  cases,  namely  for  misad justment  of  the  phase  shift 
between  the  two  voltages  feeding  the  crystals,  and  for 
misalignment  of  a  crystal . 


In  case  the  voltage  phase  shift  is  not  ^ ,  we  get  the 
following  expression  for  output  components  of  the  optical 
beam: 
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The  spectral  content  of  this  output  is  again  obtained 
through  the  "Golden  Rule"  of  Bessel  functions: 


Hence,  the  major  consequence  of  a  misad justment  in 
the  phase  shift  is  the  appearance  of  the  (-1)  order,  with 
the  following  amplitude: 


Instead  of  being  completely  suppressed,  as  is  the 
case  when  ^  '  the  residual  (-1)  order  will  be  directly 
proportional  to  the  phase  error,  and,  since  it  has  the 
same  polarization  as  the  (+1)  order  (upper  sideband) ,  it 
will  go  through  the  circular  analyzer  without  being  atten¬ 
uated.  The  harmonic  content  of  the  output  will  also  be 
made  of  the  leakage  from  the  lower  sideband  and  all  the 
(4n  +  3)  orders  which  do  not  exist  when  the  two  voltages 
are  exactly  in  quadrature. 

The  leakage  from  the  lower  sideband,  by  far  the  most 
important,  is  then: 
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Another  case  of  gj.eat  importance,  one  which  may 
degrade  very  easily  the  performance  of  the  frequency 
shifter,  is  a  misalignment  of  the  crystals.  When  the 
beam  is  not  propagating  along  the  axis  of  the  crystal, 
the  natural  birefringence  of  the  crystal  as  seen  by  the 
incoming  light  may  be  of  the  same  order  of  magnitude  or 
greater  than  the  birefringence  induced  by  the  driving 
voltage.  The  latter  is  always  less  than  a  half-wave  if 
one  does  not  want  to  damage  the  crystal,  while  natural 
birefringence  caused  by  misalignment  may  be  several  wave¬ 
lengths,  depending  on  the  thickness  of  the  crystal  and 
the  values  of  the  ordinary  and  extraordinary  indices.  In 
crystals  like  KDP,  the  difference  between  the  two  indices 
is  rather  large,  so  that  a  very  small  displacement  of  the 
axis  is  sufficient  for  a  substantial  birefringence  to 
appear.  This  problem  does  not  exist  in  cubic-type  crys¬ 
tals  like  zinc  sulfide,  which  are  isotropic  and  may  accept 
a  wide  cone  angle  of  incoming  light  without  affecting  its 
state  of  polarization. 
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In  order  to  determine  the  effect  of  crystal  misalign¬ 
ment,  we  compute  the  amount  of  natural  birefringence  ex¬ 
perienced  by  an  incoming  beam  of  light  not  propagating 
along  the  optic  axis  of  the  crystal.  This  is  done  by 
considering  the  intersection  of  the  index  ellipsoid  with 
the  plane  normal  to  the  direction  of  propagation  of  light. 
This  intersection  is  usually  an  ellipse  whose  equation 
can  be  written  with  respect  to  the  principal  indices 

nv'  nx,>  nT  crystal  and  the  direction  cosines  oi  /} 

x  y  z  '  /  '  l 

of  the  incident  light: 


After  evaluating  the  lengths  of  the  ellipse's  main 
axes,  one  finds  that  the  phase  difference  between  the  two 
linear  components  of  light  propagating  through  the  crystal 
can  be  written  as  follows: 


cos 


-  (nz  -  n,)sin9,  5 in6t 


where  and  are  the  angles  between  the  wave  normal 
and  the  two  axes  of  the  crystal.  The  additional  phase 
difference,  created  by  the  fact  that  light  is  no  longer 
propagating  along  the  optic  axis  of  tha  crystal,  introduces 
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a  set  of  harmonics  which  show  up  in  the  output  of  the 
frequency  shifter  and  degrade  its  spectral  purity.  It 
is  easy  to  see  that  when  no  voltage  is  applied  on  one  of 
the  electrooptic  crystals,  t.i\e  natural  birefringence  intro¬ 
duced  by  it  is: 

f  4>  «  13.  si  /»  .»  )a,»S 

T  >  cosr  1  v 

The  leakage  of  other  orders  than  the  hoped-for  side¬ 
band  is  directly  proportional  to  this  phase  shift,  which 
sets  a  maximum  misalignment  angle  for  both  crystals  if 
a  high  spectral  purity  is  to  be  obtained. 

Since  the  electro-optic  effect  depends  upon  the  length 
of  the  crystals,  it  is  doubled  if  the  beam  goes  twice  into 
the  crystals,  for  example,  after  reflection  upon  a  mirror. 
So,  we  have  also  tried  the  following  configuration,  with 
double  passage  of  light  through  the  crystals: 
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As  before,  the  axes  of  the  crystals  are  at  45°  with 
respect  to  each  other,  the  voltages  are  in  quadrature,  and, 
by  applying  Jones  calculus  to  the  beam  of  light  on  .it9  way 
back,  we  find  that  the  output  light  right  before  the 
circular  analyser  has  the  following  components: 

An  = 

An  *1  =  ® 

A„  -  f 


All  the  odd  harmonics  go  through  the  circular  analyser 
without  being  attenuated,  while  the  even  ones  are  eliminated 

by  it.  The  theoretical  efficiency  curves  are  plotted  on  p.30 
as  well  as  some  experimental  results.  We  can  see  that,  con¬ 
trary  to  the  preceding  single-pass  set-UD.  the  lower  side¬ 
band  is  oresent  in  the  output,  but.  for  voltages  lower  than 
the  quarter  wave  voltage,  it  is  quite  negligible,  compared 
to  the  upper  sideband,  due  to  the  linearity  of  the  first- 
order  Bessel  function  in  that  region.  It  is  also  seen  that 
the  efficiency  of  transforming  the  intensity  from  the  orig¬ 
inal  frequency  to  the  upper  sideband  is  twice  the  efficiency 
obtained  with  the  single  pass  set-up,  for  low  voltages, 
since  then: 


29 


Fig.  3.  Frequency  conversion  efficiency. 

-  -  -  =  single  pass  two  crystal  set-up 

-  =  double  pass  two  crystal  set-up 

-  =  cubic  crystal  set-up 

■**  =  experimental  data 
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4  =  K  j/t{z)  +  j,(£) 

'JcUt  2.  ^  ' 


pass, 


However,  as  far  as  the  final  efficiency  and  output 


stability  are  concerned,  it  does  not  seem  that  any  sub¬ 


stantial  gain  is  obtained  over  the  previous,  single  pass 


system.  Due  to  the  numerous  surfaces  along  the  path  of 


the  beam  (i.e.,  linear  polarizer,  quarter  wave  plate,  the 


two  crystals,  quartz  windows) ,  and  due  to  the  absorption 


by  each  of  the  components,  the  insertion  loss  is  of  the 


order  of  50%  for  one  passage  through  the  frequency  shifter. 


Hence,  the  double  pass  insertion  loss  counterbalances  the 


theoretical  gain  in  efficiency  expected  from  the  double 


pass  configuration.  Furthermore,  this  set-up  wastes  at 


least  another  half  of  the  incoming  light  by  letting  the 


beam  go  twice  through  the  beam  splitter.  Hence,  the  max¬ 


imum  efficiency  in  the  double  pass  configuration  should 


be  at  least  four  times  larger  than  in  the  single  pass  in 


order  to  achieve  the  same  power  conversion  from  the  initial 


frequency  to  the  desired  shifted  frequency.  We  have  shown 


that  in  the  best  case  (small  exciting  voltage)  it  can  be 


31 


at  most  four  times  as  large,  so  that  for  the  same  applied 
voltage,  less  flux  is  obtained  at  the  shifted  frequency 
in  the  double  pass  configuration  than  in  the  single  pass. 

Another  problem  we  have  encountered  in  the  double 
pass  system  is  the  feedback  from  reflection  upon  the  out¬ 
put  mirror  of  the  laser.  We  have  seen  previously  that  if 
the  beam  of  light  is  not  propagating  along  the  optic  axes 
of  the  crystals,  natural  birefringence  rapidly  increases 
in  value  over  induced  birefringence  and  makes  some  of  the 
light  at  the  initial  frequency  appear  in  the  output. 

If  we  want  to  keep  the  spectral  purity  of  the  output 
very  high,  for  example,  with  the  energy  being  in  the  first 
order  by  more  than  a  certain  percentage,  say  99.29%,  so 
that  the  amplitude  in  the  zeroth  order  is  less  than  1%  of 
the  amplitude  of  the  first  order,  we  have  shown  above  that: 


^  0.0/ 


where  Q  is  the  thickness  of  the  crystals  and  the 

maximum  misalignment  between  the  crystals'  axes  and  the 
beam.  With  the  actual  characteristics  of  the  set-up,  we 
find  that: 

=  3*  ,0~l  rdc'ien 
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Therefore,  the  beam  has  to  come  back  as  much  as  possible 
alonq  the  same  way  it  came  in.  Care  must  be  taken  that 
the  divergence  of  the  beam  from  the  laser  does  not  exceed 
this  limit,  which  is  the  case  for  the  laser  we  are  using. 
If  the  divergence  of  the  laser  beam  was  larger,  one  should 
match  the  beam  to  the  crystals  by  putting  the  beam  waist 
on  the  mirror  as  close  as  possible  to  the  crystals. 

When  this  alignment  is  done,  much  fluctuation  can 
be  noticed  in  the  output  of  the  frequency  shifter,  which 
can  be  explained  by  the  fact  that  the  mirror  of  the  laser, 
the  beam  splitter  and  the  final  mirror  act  like  a  low  re¬ 
flectance  Fabry-Perot,  internally  modulated  in  frequency 
by  the  electro-optic  crystals: 


mt77T> 


rt 
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As  seen  previously,  the  frequency  response  of  the 
double  pass  frequency  shifter  (without  any  feedback  from 
the  laser)  to  a  quasi-monochromatic  beam  of  light  can  be 
written  as  follows: 


Ll t 


Z  A 


t(u.  *  "HJ 


where  the 


coefficents  are  given  on  page  29.  At  each 


passage  back  through  the  system,  each  harmonic  gives  rise 
to  an  infinite  series  of  outputs,  and  we  find  that  the 
output  beam  after  the  beam  splitter  is  made  up  of  harmonics 
with  the  following  amplitudes: 

K~tZ  pu,'E£z...  4  4  4 ... .  "r 

k  n-o  r  7  r  r  !  r 

where:  f)  ^  -  £ mpji /</</*  of  fhc  harmonic. 


tt  ^  s.  plihelc  bransmi  Mancc.  and  rtfltctsnc*. 

of  flit.  sp/i fft-r. 

if  =  ^5  Zl 

I  > 

p  +  =  k 

For  voltages  lower  than  the  quarter  wave  voltage  of 
the  crystals,  the  only  non-negligible  output  for  the  double 
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pass  frequency  shifter  without  feedback  is  the  first 
order,  so  that  the  amplitude  of  the  optical  field  at  the 
outset  of  the  beam  splitter  is,  taking  the  feedback  into 


account : 


Any  change  in  the  distance  [  between  the  laser  and 
the  mirror  M  will  modify  the  length  of  the  cavity,  hence 
,  and  modulate  both  the  amplitude  and  phase  of  the  out¬ 
put  optical  field  represented  by  the  phasor  c/^r.  In  order 
to  eliminate  this  problem,  we  have  had  to  put  an  optical 
isolator  (  made  up  of  a  quarter  wave  plate  followed  by  a 
linear  polarizer) ,  which  prevents  the  light  reflected  from 
the  output  mirror  of  the  laser  from  being  re-injected  into 
the  frequency  shifter. 


III.  A  METHOD  FOR  OPTICALLY  CORRECTING  CRYSTAL  AND 
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GLASS  LASER  RODS 

1 .  Review  of  Earlier  Progress 

As  described  in  earlier  contract  reports,  the  sequence 
of  steps  in  manufacturing  an  optical  correction  plate  is 
as  follows: 

(1)  The  distortion  to  be  corrected  is  determined  by 
obtaining  an  interf erogram  of  a  wavefront  which 
has  been  distorted  by  passage  through  the  laser 
rod  in  question.  This  interferogram  is  essential 
ly  a  contour  map  of  the  distorted  wavefront,  with 
contours  spaced  by  a  half  wavelength,  and  also 
represents  a  contour  map  of  the  optical  thickness 
of  the  desired  correction  plate  (with  high  spots 
switched  to  low  spots,  etc). 

(2)  The  correction  plate  is  made  by  grinding  a  glass 
plate  to  match  the  contour  map  defined  above, 
except  that  a  scale  factor  M  will  be  included; 
e.g.  ten  waves  of  distortion  may  be  scaled  up  to 
an  appreciable  fraction  of  a  millimeter  so  that 
the  contouring  can  be  accurately  accomplished 
using  an  optical  milling  machine. 
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(3)  Finally,  the  "magnified"  correction  plate  is 
scaled  back  down  to  the  proper  magnitude  by 
immersing  the  plate  in  a  fluid  whose  refractive 
index  very  nearly  matches  that  of  the  glass 
correction  plate.  An  index  mismatch  of  An 
produces  a  scale  factor  M  =  ]^,n  .  This  tech¬ 
nique  of  using  near  index -matching  provides  a 
number  of  additional  advantages,  the  most 
important  of  which  is  that  it  does  away  with 
the  necessity  for  optically  polishing  the 
correction  plate  since  the  surface  roughness 
of  the  ground  correction  plate  is  scaled  down 
by  the  same  factor  M  .  This  prediction  was 
confirmed  experimentally  at  an  early  stage  in 
this  work. 

To  test  the  validity  of  this  entire  concept,  we  decided 
to  pick  a  piece  of  window  glass,  exhibiting  about  20  fringes 
of  random  deformation,  and  make  a  correction  plate  for  it 
in  the  manner  described  above.  The  sample  to  be  corrected 
was  approximately  2  inches  square.  The  interferogram  was 
photographed  using  a  Twyman -Green  interferometer  with  a 
546  iR  mercury  light  source.  This  photographic  interferogram 
was  then  digitized  using  an  HP  digitizer  as  described  in 
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detail  in  earlier  reports.  The  digitized  data  were  then 
used  to  generate  a  paper  tape  which  in  turn  was  used  to 
control  an  optical  milling  machine.  This  milling  machine, 
which  was  simply  a  modified  metal-working  machine  owned 
by  the  Frank  Cooke  Optical  Company,  could  be  controlled 
in  cutting  depth  during  constant  feed  by  a  paper  tape 
reader.  The  transfer  of  data  from  the  original  digital 
form  to  the  proper  paper  tape  format  proved  to  be  a  very 
lengthy  and  cumbersome  procedure  -  but  would  of  course  be 
much  simpler  the  next  time  around.  The  remainder  of  this 
report  describes  the  operation  of  the  tape -controlled  mill¬ 
ing  machine  in  some  detail,  and  assesses  the  quality  of 
the  correction  plate  which  was  made. 

2 .  Production  of  Scaled  Correction  Plates  using  a  Tape- 

Controlled  Milling  Machine 

A.  Machine  control 

As  described  earlier,  by  obtaining  a  digitized 
two-dimensional  mapping  of  the  interf erogram  we  provide 
ourselves  with  a  digital  description  of  the  distorted 
wavefront,  and  hence  of  the  desired  correction  plate.  The 
index  matching  scale  factor  M  is  then  applied  to  generate 
the  surface  coordinates  of  the  machined  correction  plate. 

In  addition,  in  order  to  overcome  the  problem  of  back-lash 
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In  the  milling  machine,  the  data  were  transferred  to  a 
tilted  coordinate  system  (Fig.  4)  so  that  the  motion  of 
the  cutting  tool  was  always  in  the  same  direction. 


MOTION  OF 


Fig.  4  Set-up  for  machining  correction  plate 

Figure  5  shows  the  relative  positions  of  the  cutting 
tool  and  the  glass  plate  at  the  beginning  of  a  cutting 
"scan'1. 
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DESIRED  SURFACE 
CONTOUR 


PATH  OF  POINT  Q 
RELATIVE  TO 
CORRECTION  PLATE 


POINT  OF 
CONTACT 

Fig.  5  Relative  positions  of  cutting  tool  and  glass  plate 
at  the  beginning  of  a  cutting  "scan” 


^s  shown,  the  contour  of  the  machined  plate  is  actually 
letermined  by  the  path  of  the  point  Q  ,  which  is  the  point 
3 f  contact  between  the  tool  and  the  glass  plate.  At  the  end 
3 f  a  scan,  the  tool  is  returned  to  the  starting  position 
and  the  correction  plate  is  manually  stepped  in  a  direction 
perpendicular  to  the  plane  of  Fig.  5.  Control  of  the  machine 
is  then  returned  to  the  tape  reader  and  the  ne^t  adjacent 
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cutting  scan  is  started.  In  producing  the  paper  machine 
control  tapes  we  start  by  storing  both  the  scaled-up  and 
transformed  data  sets,  and  also  some  control  signals  :e- 
quired  for  the  milling  machine  operation.  This  information 
was  stored  in  the  disc  of  an  IBM  360,  and  then  transferred 
at  a  teletype  terminal  to  paper  tape.  Since  the  code  for¬ 
mat  of  the  teletype  is  the  ASCII,  we  had  to  subsequently 
convert  it  to  the  EIA  format  used  by  the  milling  machine. 

This  was  done  with  the  aid  of  an  HP  Computer  Data  Acquisition 
System  --  HP-2748A/Tape  Reader,  HP-2440A  Analogue -Digital 
Interface,  and  HP-7900A  Disc  Drive.  As  mentioned  earlier, 
this  was  an  extremely  awkward  arrangement  and  could  certain¬ 
ly  be  circumvented  in  the  future. 

B.  Design  of  the  cutting  wheel. 

Although  diamond -bonded  glass -cutting  tools  can 
now  be  made  in  virtually  any  shape  and  with  any  desired 
size  diamond  grit,  we  ended  up  using  one  that  Frank  Cooke 
had  on  hand.  Following  are  some  of  the  considerations  in 
designing  the  optimum  shape  tool.  The  tool  is  a  section 
of  a  sphere  as  shown  in  Fig.  6. 
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From  the  point  of  view  of  wear  on  the  cutting  tool,  which 
is  far  more  subject  to  wear  than  a  conventional  metal¬ 
cutting  tool,  the  cutting  surface  area  -  and  hence  the  tool 
radius  R  -  should  be  as  large  as  possible.  On  the  other 
hand,  fine  control  over  the  surface  figure  of  the  correction 
plate  requires  a  small  radius  tool.  Fig.  7  shows  a  number 
of  parameters  of  the  problem.  In  particular,  the  "residue" 
is  defined  as  the  height  of  the  cusp-shaped  residual  material 
on  the  correction  plate  between  cutting  scans  spaced  by  a 
distance  h  .  This  residue  must,  of  course,  be  small  enough 
to  be  reduced  well  below  a  wavelength  after  index  matching. 
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It  can  easily  be  shown  that  the  radius  of  the  tool  must 
satisfy  the  following  inequality  in  order  to  leave  a 
residue  no  larger  than  d  ,  with  a  spacing  between  cutting 
scans  of  h  . 

R  >  h2/8d 


RESIDUE  d  LEFT 
BETWEEN  CUTTING 
SCANS 


Fig.  7  Machining  of  the  correction  plate 


The  thickness,  T  ,  of  the  cutting  wheel  is  not  critical 
but  should  certainly  be  large  enough  to  ensure  that  the 
tool  will  not  vibrate  or  change  shape  during  grinding. 
Obviously,  the  tool  thickness  must  be  greater  than  the 
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spacing  h  between  cutting  scans. 

Strictly  speaking,  the  surface  figure  produced  by 
cutting  at  the  point  of  tangency  between  the  tool  and 
the  correction  plate  is  not  identical  with  the  locus  of 
the  cutting  tool  axis  relative  to  the  plate.  This  dis¬ 
crepancy  could  have  been  taken  into  account  in  producing 
the  control  tape,  but  was  in  fact  small  enough  to  safely 
neglect;  the  smaller  the  radius  of  the  cutting  tool,  the 
smaller  the  discrepancy. 

3 .  Assessment  of  Finished  Correction  Plate 

The  performance  of  the  completed  correction  plate  is 
illustrated  in  Figs.  8-11.  A  number  of  features  are  im¬ 
mediately  apparent: 

1.  Because  (i)  we  were  unable  to  obtain  the  proper 
index  matching  fluid,  and  (ii)  the  residue,  d  , 

(see  Fig.  7)  was  somewhat  excessive  due  to  the 
fact  that  we  spaced  the  cutting  scans  too  far 
apart  in  order  to  keep  the  machining  costs  down 
to  a  reasonable  figure,  the  interference  fringes 
in  the  index-matched  correction  plate  have  about 
a  half -wavelength  of  deformation  due  to  the  residue. 
(This  is  why  the  fringes  in  Figs.  9-11  are  "jagged".) 
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Fig.  8  Distorted  wavefront  to  be  corrected  (2"x2") ; 

approximately  20  fringes  (Twyman-Green  interferogram) 


^^•8*  9  Interferogram  of  imperfectly  index -matched  correction 
plate 
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Fig.  10  Interferogram  of  corrected  wavefront  (approximately 
4  fringes) 


Fig.  11  Interferogram  of  imperfectly  index -matched 

correction  plate  with  reference  wavefront  slightly 
tilted. 
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Fig.  12  Photo  of  one  edge  of  correction  plate,  showing 
grooves  and  residue  left  by  machining  operation 
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2.  The  correction  plate  was  successful:  when  placed 
in  series  with  the  plate  to  be  corrected,  a  mark¬ 
ed  improvement  in  the  quality  of  the  transmitted 
wavefront  was  obtained.  The  uncorrected  wave- 
front  had  about  20  fringes  of  deformation,  while 
the  corrected  wavefront  had  just  a  couple  of 
fringes.  The  improvement  is  particularly  evident 
in  Fig.  11,  where  some  tilt  was  introduced  in  the 
interferometer  to  more  clearly  display  the  nature 
of  the  residual  wavefront  deformation. 

We  feel  that  this  work,  which  consisted  of  (i)  measur¬ 
ing  the  deformation  to  be  corrected,  (ii)  digitizing  the 
deformation  information,  (iii)  preparing  a  paper  tape  to 
control  an  optical  milling  machine,  (iv)  machining  a  scaled- 
up  correction  plate,  and  (v)  using  the  correction  plate 
with  the  appropriate  index -matching  fluid  to  attempt  to 
correct  the  original  deformation,  has  been  successful.  We 
have  clearly  demonstrated  the  feasibility  of  this  approach 
to  correcting  arbitrary  wavefront  deformations.  With  our 
current  capability  we  could  easily  correct  wavefront  dis¬ 
tortions  or  up  to  50  wavelengths,  to  less  than  a  wavelength. 
Our  approach  to  this  problem  would  be  particularly  well-suited 
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to  the  correction  of  large  numbers  of  optical  components, 
since  initial  programming,  digitization,  and  optical 

machining  can  all  be  done  relatively  cheaply  and  routinely 
after  the  initial  set-up. 


APPENDIX 


Program  for  Making  a  Correction  Plate 
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/,iatfor 

c  PROGRAM  KOw  mAnIng  A  CORRECTION  PLATE 

C  CMkN'  INSTITUTE  OF  TPTlCSr  UNIVERSITY  OF  ROCHESTER 

DIMENSION  A(«0,RO),Yl«0),XX(30),YY(10l>,XP(101),ZP(101,101) 
®FAl.*n  A  r  Y  ,  DET 

l  OGICALM  TAH/210/,MlNUS/'-»/ 

rfao,nl,ll»ox,dy 

C  nl*NumuER  OF  St  A  n  LlNFS  ACROSS  InTERFEROGRAM 

C  LI  sNUMHKR  OF  INTERPOLATED  POINTS  ALONG  EACH  SUAN  LINE 

C  DXMNTfRVAL  OF  INTERPOLATION 

C  0V=INTERVAL  IN  Y-DIRECT1UN  BETWEEN  SCAN  LINES 

PRINT  101 

101  FORwa  r(  l  ix,  'NL. ',  11X,  'LL'  »BX,  '0X»,16X,  'DV»  ) 

PRINT, NL,LL, OX, DY 

00  102  L*1,LL 

102  xP(L)sDX«(L-l ) 


BEGIN  MAPPING  IN  X-DIRECTION 
DO  100  Ms 1 , NL 

REAO.Nfj,  (XX(K),KS]  ,NO),  (YV(K),Ksl,ND) 

NOsNijmPFR  OK  ORDERED  DATA  POINTS  ALONG  EACH  SCAN  LINE 
XY(H)*X-COORDINATE  OF  kTh  data  point 
vv(< )sZ-cnjRDlNATE  OF  KTM  DATA  POINT 
on  1  LSI, NO 
1  XY(L)sO,!j*XXCI_} 

T°  Pt*L  VALl'E  UF  XX*Ll  DUE  TO  ENLARGEMENT  OF  INTERFEROGRAM 

On  ?  K s | f  MO 
?  y c l isY(1)ayy(kj 

NsND 

IF(N,GT,«0)  NS«0 
on  5  I =? , n 

Y ( 1 ) SO , 
on  s  k»j ,no 
IF(T-Jb)  4,4,3 
3  I F  (  xx  (O  .LT.0.01  )  GO  TO  5 
«  X(I)sY(T)*xXfK)**(l-i)*YY(R) 

5  continue 

A ( 1 , 1 )SN0 
Do  13  J»2,N 
All,  J)sO, 

A  C  J , N ) SO  t 
GO  11  K  s  i  #  Nf) 

IF(J-«S)  7,7,6 

6  IF(XX(k),lT,P,1)  GO  TO  8 

7  A  Cl ,J)sA(] ,J)  +  XX(KJ**(J«1) 

P  JF(J*N-36)  10,10,9 

9  IF(xx(K).LT,o.d  GO  TO  11 

10  A(J,nJ«A( J,n)+XX(R) a*cj»n»2) 

11  CONTINUE 
JJ.J-1 

DO  12  Isl,JJ 

LlstPl 

LJsJ-I 

12  A(LI,LJ)3A(1,J) 

IF(J.EU.N)  GO  TO  14 
JJBN-J 
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DP  H  l=i,jj 
K  I  s  J  ♦  I 

KJSM-I 

13  A(KT,KJ1sA(J,N) 

1«  CALL  f>MATE<l(A,N,Y,OF.T) 

C  'in*  V(n=CnPrnciENT6  OF  POLYNOMIAL  (LOW  TO  HIGH) 

ZP(m, 1 )sy(i J 
DO  15  L=2,LL 
7P(m,L)=Y(N) 
no  15  1=2, n 

15  7P(*»I  )=ZP(M,L)*XP(L)*Y(N  +  l-n 
100  CONTINUE 

L>n  ?01  K  =  1,NL 
20  1  X  X  (  K  )  =1)  Y  *  (  K  “  1  ) 

N  =  NL 

IF(N,GT,«0)  N  =  «0 
C 

C  BEGIN  "APPING  In  Y-DI«ECTI0N 

L'O  200  M= 1 , LL 

Y(1)=0. 

On  21  K=1,NL 

21  Y  ( 1  )  =  Y(U*ZP(K,m> 
no  pa  1=2, N 
v(I)=0. 

OM  ?U  K=1,NL 
IMT-35)  21,23,22 

22  IF(xxU).l.T.O.Ol)  GO  TO  2« 

23  Y(I)=V(T)fXXCK)**(I-l)*ZP(K,M) 

2u  CONTINUE 

A (1 , ll=NL 
no  3?  J=2,N 
A(1,JJ=0, 

A (J,N)=0, 

On  30  *<  =  1,NL 
JF(J-3S)  26,26, 25 
2S  IF(XX(k),LT,0,1)  GO  TO  27 
2h  A(l,Jl=A(l,J)fXX(K)AA(J-l) 

27  IP  (  '*N-3b) 

2«  IF(xx(K).|,T.0.1  )  GO  TP  30 

29  A  (  J,M=A(  J,U)«xX(K)**(J*N-2) 

30  CONTINUE 
JJsJ-1 

no  3i  i=i, jj 
11*1*1 
LJsJ-I 

3 1  A(LI,LJ)=A(1,J) 

T  F  ( J ,  E  fg ,  n  )  Gn  TO  33 

JJSN-J 

nn  32  1=1 , J J 
KT=J*I 

* JSN-I 

32  A(Kl,«J)=A(.?,N) 

33  CALL  0matFn(A,N,Y,UET) 

C  YdlsCOFFFKIENTS  OF  POLYNOMIAL  (L0w  TO  HIGH) 

C  *  A  v  F  l  t NG  Ths  5061  A  ,  SCALING  FACTOR  M  =  200 

C  lo  In  THE  UNIT  OF  1/10000  INCH 

ZP(1 ,m)=Y(1 )*2i,5 
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Dr  35  L*2,LL 
7f*lL  ,*)sY(N) 

Dn  34*  1*2* v 

ZP(L,«)*lPCI.,P)»xP(L)*VCN*|*n 

3S  ZP(L.^)*ZPCL*M)*2itS 
t  ?p  Im  T*fc  UNIT  nt  1/10000  Inch 

200  CPNTINUF 
C 

bOp  FnkMAT(Il,Al,I3) 

SOI  FnkMAT(J2,Al,I3) 

SO?  f  C'hN*T(!3,*i,i3) 

S03  ^0K»-4TCI<,,a1,13) 

600  F^*-«AT(I1  ,41  **l,Ifl) 
bOi  FnNf«AT(J2, 41 ,41,14) 

602  F  PF  *«  a  T (13*41*41*14) 

603  Format  (la, 4]  ,4i,i<ij 
SOs  Fnr-^AT(Tl*Ai,Ii*«|*I3) 
b0«  Fn6M*TiI<i,Ai,i2,*i,i3) 

SOQ  Fr6*AT(I<i,Al,Il,A!,I3) 

510  F  OF  m  a  T (T3*a1 *  12  *  a  1 *13) 
bll  Fnt.-ATC!3,A1,II,A1*I3) 

512  Fn**»'*T(I2,ti,I2,Al*I3) 

513  Fnfc^AT (I2*41*I1*A1»I3) 

5M  FOf>**AT(Tl  *41*12,41*13) 

515  FOMATH  1,  A  1,  A  1,41,15,41,41  #A|,  12) 

516  Fps*'aT(T2,AI,a1,aI,I5,A1,41,41,I21 

517  FOM-aT (I  3, A  1,4 l, A  1,15, A  1,41, 41, 12) 

51  e  F  OP- AT  U«,A1,  A  1,41,15,  A  1,41, 41, 12) 

520  Fn„-AT(Il,Al,Al,J<|,Al,Al,Al,»,,Il,li,Iif) 

521  F  P664V  (12, 41  ,  A  1,I<I,  A  1,41,41,41,11,11,12) 

522  FOh^AT (13,4!, 41 ,!«,4 1,41 ,41,41,11,11,12) 
5?5  FPt*-4T(I«,4l,Al,Ill,Ai,41,41,Al,Il,ll,I2) 

*-nx*60o 

NHYS5100 

C  NMX  SOY  FOR  TRANSLATION 

S380. 0523a 
C3*0.QRP63 

C  S3  C3  for  TRANSFORMATION 

I  h*S5 
IS*0 
I  T=6 

C  IH  IS  IT  FOR  CONTROL  SIGNALS 

NO*0 

C  NO  is  T WE  SF«UtNCE  NUMBER 

DO  4QR  1*1,48 

c 

C  TRANSFORMATION  and  TRANSLATION 

rn  406  J*i ,ll 

VV(J)*Kr.VtKP(J)*iooOO,AC34ZP(I,J)tS3 

lr*VY(J) 

IF( YV(J)-lY-0,5)  401,401,402 

401  VV( J)*lY 
GO  TO  403 

402  YV(J)«lYtl 

4  03  ZP C I * J )*NO*.( xP ( j  j *  1 0000 * *(«S3)  +  ZP( I, J ) *C3) 
I?*ZP(I,J) 
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ifizpu, j)-iz-o,s)  404,404,405 

aOa  ZP(I,J)sI7 
GO  TO  <106 
«0b  /P(I,J)sIZM 
406  CONTINUE 

nosnoM 

ixs/p(i,n 

If-  (N0,GT,999)  UO  TU  <113 
IF(NU,GT,99)  GU  IP 
IF(Mi.GT,9)  GO  TO  ail 
1 Tt ( 10,SliO)  NO, TAB* IX 
go  to  at  a 

ail  *9ITFClO,501>  MO, TAB, IX 
Go  To  aia 

<*1?  MJlTF  (  10, So?)  NO,TAB,IX 
GO  Til  414 

413  fcOlTt (10,503)  Ml, TAB, IX 
ala  I'OsMHi 
TVsVY(l) 

IF(MU,GT.999)GU  TO  417 
IE(Ml,GT,99)Gfl  TU  416 
IE  (M).GT  ,9 ) GO  10  415 
,.«lTE(10,600)Nu,Hh,TAB,IY 
GO  TO  aia 

alS  a91TE(  10,601 ) no, TAB, TAB* I Y 
Gn  To  al9 

ai6  *PmCi0,60?)M),TAB,TAB,lY 
GO  TO  a  1 9 

a  1  7  *91  U(  i  0,603)NIJ,TAH,TAB,IY 
419  no  4.30  Js?,LL 
r  nswi+i 

IXSZP(I,J)-ZP(I,J-1) 

IVSYYf J)-YY(J-1 ) 

IE (N U, GT, 999. AND. IX.GT.0)  GO  TO  423 
IE (k ".GT .999.A1 o, iX.tfl.O)  GO  TO  424 
IE(M.,GT.99,ANU,lX,GT,0)  GO  TO  425 
IE(*0,bT,99,AM»,IX,FUtD)  GO  TO  426 
IF  (M..GT  .9,  AvO.  IX.GT.  0)  GO  To  927 
IE  CNO.GT.a,  A\r.,lx.t0,0)  GO  TO  428 

IE(IX.GT.O)  GO  T(J  429 

»'BlTl  ( 1  0,505)  NO, TAB, IX, TAB, IV 

GO  TO  430 

4?3  *B1TE(10,S0A)  no,taH,IX,TAB,IY 
GO  TO  430 

42a  *»ITr ( 10, 509)  NO, TAB , I x , T AR, I y 

Go  TO  430 

425  * BITE ( 10,510)  fn,TAB,Ix,TAH,lv 

GO  TO  430 

'< 2 6  *»ITt  (10,51  1  )  NO, TAR,  IX, TAB,  IV 
GO  TO  430 

42  f  *  °I TF { l 0,512)  Ml, TAB, IX, TAR, IY 
GO  To  430 

42fl  *  R  1  Tt  (10,513)  NO, TAB, IX, TAB, IY 
GO  Tn  430 

4?9  *«lTtf 10,514)  NO, TAR, IX, TAB, IY 
a30  CONTINUE 
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NO=NOM 

IV5VY(LU 

IF(N0,GT,999)  GO  TO  993 
IF  (Ml),  GT  ,  99  )  GfJ  TO  992 
IF(MJ,GT,9)  GO  TU  99l 

*RlTF<10,5l5)  NO,  TAB#  TAB, MINUS , I Y , TAB# TAB, TAR,  IH 
GO  TO  950 

991  fcPlTt(10.5l6)  NO,TAB,TAB,MINL,S,lY,TAH,TAB,TAR,IH 
Gn  TO  950 

9«2  wrITE(10,517)  NO, TAB, TAB, MINUS, IT, TAB, TAB, TAB, IH 
GO  TO  950 

993  *RlTE ( 10,5t«)  NO,  TAB,  TAB, MINUS,  IY#TAB, TAB, TAB#  IH 

950  NOSNOfl 
IX=ZP(I,LL) 

IFUtEu,«a)Ils2 

I F  C  NO,GT  f  999 }  GO  TU  953 
IF(N0,GT#99)  GU  Tf)  9b? 

IFCNU.GT.9)  GM  TU  951 

"RITKl 0,520)  NO, TAB, MINUS, I X, TAB, TAB, TAB, TAB, IS, IT, IH 
GO  TU  960 

951  WR I  TF ( 1 0 , 52 1 )  NO, T AB , MInUS , I X , T AB , T AB , T AB , TAB , I S, I T , IH 
GO  TO  960 

95?  WRITE (10, 52?)  NO, TAB, MINUS, IX, TAB, TAB, TAB, TAB, IS, IT, IH 
GO  TO  960 

953  WRITF(10,523)  NO, T AB # M INUS , I X , T AB » T AB » T AB# T AB# IS, I T , I H 
960  CONTINUE 
999  CONTINUE 
STOP 
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